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ABSTRACT 

The gas from which stars form is magnetized, and strong magnetic fields can efficiently transport 
angular momentum. Most theoretical models of this phenomenon find that it should prevent formation 
of large (> 100 AU), rotationally-supported disks around most protostars, even when non-ideal 
magnetohydrodynamic (MHD) effects that allow the field and gas to decouple are taken into account. 
Using recent observations of magnetic field strengths and orientations in protostellar cores, we show 
that this conclusion is incorrect. The distribution of magnetic field strengths is very broad, and 
alignments between fields and angular momentum vectors within protostellar cores are essentially 
random. By combining the field strength and misalignment data with MHD simulations showing that 
disk formation is expected for both weak and misaligned fields, we show that these observations imply 
that we should expect disk fractions of ~ 10 - 50% even when protostars are still deeply embedded in 
their parent cores, and even if the gas is governed by ideal MHD. 

Subject headings: magnetic fields — protoplanetary disks — stars: formation — stars: protostars — 
stars: winds, outflows 



1. INTRODUCTION 

Stars form via the collapse of dense, magnetized cores 
of interstellar gas. As the gas contracts, it spins up to 
conserve angular momentum and twists magnetic field 
lines. This launches torsional Alfven waves that transport 
angular momentum away from the densest, collapsing 
region. In cores that have aligned field and rotation axes, 
and flux-to-mass ratios within a factor of a few of the 
critical value (defined as the value at which the magnetic 
field is able to prevent collapse entirely), this magnetic 
braking is so effective that Keplerian disks do not form 



(Galli et al. 2006; Price &: Bate 2007; Mellon fc Li 2008 



'. iennebelle fc (Jiardi|2009| |Kr asnopol sky et al.|2012 j. For 
fields that are stronger (but that are still weak enough to 
render the core magnetically supercritical), the outcome 
is no disk at all; for weaker fields, the result is a sub- 
Keplerian "pseudodisk" supported by magnetic pressure 
rather than rotation. 

Since Keplerian disks and the planetary systems they 
produce are observed to be ubiquitous aroun d optically- 



visible young stars (e.g., Haisch et al.||200l" ), this result 



presents an obvious problem. A number of possible solu- 
tions have been proposed. One is that disks do not appear 
until the majority of the surrounding protostellar core has 
been accreted, at which point the inertia of the envelope 
into which the twisted fields deposit angular momentum 



Mellon & Li 


2009 


Krasnopolsky et al. 


2012 



In this 



when the envelope is cleared and protostars become opti- 
cally visible. This proposal faces two severe challenges, 
however. First, observations now directly demonstrate 
that extended disks, at least some of them Keplerian, 
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large envelopes ( 


J0rgensen et al.||2009[ |Enoch et al.||2009 
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And second, 



the magnetic field lines threading the pseudo-disk still 
connect to the much larger surrounding molecular cloud, 
and it is unclear why magnetic braking would cease to 
operate. 

Another possible solution is that non-ideal MHD effects 
allow the magnetic field to decouple from the gas, pre- 
venting efficient magnetic braking. A number of authors 
have investigated whether ion-neutral drift (ft 



[2009| 



IDufiin & Pudritz 2009: Li "eTaLlfeOl 



M ellon fc Li 
or (Jhmic 



issipation (Dapp fe; Basu||2010| |Machida et al.||2011 Li 



et al.|20il I might allow Keplerian disks to form, and have 



found that either they fail to produce disks at all, or they 
produce disks that are only ~10 AU in size. In contrast, 
the observed disks even around deeply embedded proto- 
stars are an order of magnitude larger (|j0rgensen et al.| 
20091 [Enoch et al.||2009l |Takakuwa et al.| |2012| |Tobin| 



et al. 2012p . Turbulent reconnection might also allow 
disk formation (Santos-Lima et al. 2012), but thus far 
this has been demonstrated to be effective only in the 
presen ce of supersonic turbulence. In contrast, observa- 
tions (|Kirk et al.||2007[ " 
l[2oW 




Andre et al 
Pineda et al. 112010) and t 



2007 



tieory 



Rosolowskyl 
( juftner et aT| 



show that typical low-mass stellar cores are at most 



trans-sonically turbulent. 

A third option, which we investigate here, is that disk 
formation might be enabled by a combination of low mag- 
netic field strengths and misalignment between the fields 
and the angular momentum vector of the protostellar 
cores they thread. Simulations of misaligned cores indi- 
cate that extended Keplerian disks can form if the field 



4 Note that |Maury et ah] | |2010| have searched for fragmentation 
or other circumstellar structures around class sources on few 
hundred AU scales, and report a null result. However, due to their 
interferometer settings, they are insensitive to disks with masses of 
~ 0.1 Mq or less. 
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is weak enough and the angle between the field and the 
rotation axis is large enoug h dHennebelle fc Ciardi | [2009 



Ciardi k Hennebelle||2010| jjooset al.||2012j . While this 



possibility is promising, up until now it has not been 
possible to determine if this is a viable option for disk 
formation in general, because the distribution of core 
magnetic field strengths and field/rotation misalignments 
was unknown. However, recent observations have changed 
this situation, providing good statistical measures of both 
quantities. 

In the following sections we first discuss the observa- 
tional evidence regarding magnetic field strengths, and 
then that regarding field/rotation alignment. We then 
combine these observations with simulations in order to 
calculate disk fractions. Finally, we discuss the implica- 
tions of our results. 

2. OBSERVATIONAL INPUTS 
2.1. Magnetic Field Strengths in Protostellar Cores 
Field strengths can be measured via the Zeeman effect 



( Crutcher etaL|2010[ ), although with current sensitivities 
it is only possible to measure the line-of-sight component 
i?i os °f the magnetic vector B. Hence, only lower limits 
on the total magnetic field strength -B to tai, and therefore 
lower limits on the magnetic flux-to-mass ratio <&/M, 
may be measured directly. If gravity is balanced by 
magnetic support, $/Af is said to be critical. Hereafter 
all $/M values are normalized by the critical value, so if 
$/M < 1 for a cloud it is supercritical and the magnetic 
field is too weak to support the cloud against gravity. For 
a sufficiently large sample of clouds one would expect 
some B to be approximately along the line of sight, so if 
subcritical clouds existed some observed $/M lower limits 
would be above the critical value. However, evaluation 
of availab le surveys of Ze eman observations of OH and 
CN lines (Crutcher 2012) has shown that there are no 
directly measured subcritical <&/M in molecular clouds; 
that is, for all Zeeman measurements of B\ OSl Q/M < 1. 

For a sample of Zeeman observations that includes 
both detections and non-detections, one obtains the PDF 
(probability density function) of B\ osl P(B\ os ). There are 
two factors that determine P(Bi os ). One is the distribu- 
tion of angles between the direction of B and the line 
of sight along which -Bi os is measured. For a sufficiently 
large sample of clouds, one assumes that the angles are 
randomly distributed. The other factor is -P(-Btotai)- Re- 
gardless of the functional form of P(-B to tai), the observed 
distribution P(Bi os ) implies that the mean value of -Btotai 
in a s ample is approximately 2 t imes the mean value of the 
-Bi os (Heiles & Crutcher 2005). (A delta- function form 
for P{B tota i) gives (.Btotai) = 2.0 x (|-Bi os |), and a flat 
-P(-Btotai) gives (-Btotai) = 1-9 x (I-Bosl).) Hence, the most 
straightforward analysis of a set of Zeeman measurements 
is to compute the magnitude of the mean of the -B to tai as 
(-Btotai) = 2 x (|-Bios|). Note that -Bi os is signed, with the 
sign indicating whether B points toward or away from 
the observer. 

A more sophisticated Bayesian analysis of the observed 
distribution of P(Bi os ) further indicates that the most 
probable fu nctional form for P(-B to t a i) is flat (Crutcher 
et al.|2010[ ); that is, at any density hh a sample ol molec- 
ular clouds would have -B to t a i ranging with equal proba- 
bility from near zero up to a maximum -B max , with that 



maximum scaling as a power law in nu- Crutcher (2012 1 
reviewed studies of ($/M) in a large sample ol molecular 
clouds and concluded that the mean flux-to-mass ratio 
($/M) sa 0.5. Combined with the results that -P(-Btotai) 
is flat, this implies that the most likely true distribution of 
$/M values is flat, with a minimum of and a maximum 
of $/M = 1. 

2.2. Field/Rotation Misalignment 

The distribution of magnetic field strengths provides 
one of the two observational inputs required to calculate 
the disk fraction. The other is the distribution of rel- 
ative angles between core magnetic fields and angular 
momentum vectors. The direction of a core's angular 
momentum vector can be inferred from the orientation of 
its outflow, since outflows are launched perpendicular to 
circumstellar disks. Field directions in the plane of the 
sky can be measured via dust polarization. Under almost 
all circumstances, spinning dust grains are expected to 
align themselves with their long axes perpendicular to 
the magnetic field ( Lazarian||2007 1 , leading to both polar- 
ized emission and absorption. In dense cores, which are 
opaque to background starlight, mapping the polarized 
thermal emission from magnetically aligned dust grains is 
the most practical means of studying the field morphology 
on small scales. 

Large-scale sub-mm dust polarization maps with reso- 
lutions of ^20" have been obtained toward many sources, 
but significantly higher angular resolution is required to 
study the field geometry in the densest regions of pro- 
tostellar cores; this requires interferometric polarization 
maps. Until recently only about a dozen sources had 
been mapped (e.g. NGC 1333-IRA S 4A: |Girart et al~ 
[20061 IRAS 16293: |Rao et al.|[2009l and Orion KL: |Rao| 
|et al.|1998| Tang et al.||2010[ ), which was too few to allow 
statistical analysis. -However, the 1 mm dual-polarization 
receiver system at CARMA (the Combined Array for Re- 
search in Millimeter-wave Astronomy) has enabled high- 
resolution ma pping of dust po larization toward many 
more sources. Hull et al. (2012) present results from the 
TADPOL survejQ a CARMA key project. They focus 
on nearby, low-mass protostars that have clear outflow 
and magnetic field orientations. 

The TADPOL results strongly rule out a scenario where 
magnetic fields and disks are aligned. Rather, their results 
are consistent with a random distribution of orientations 
between the field and the outflow. The data are even 
consistent with a preferential anti-alignment between out- 
flows and magnetic fields. If one assumes that outflows 
emerge along the rotation axes of circumstellar disks, as 
is expected in magnetocen t rifugally launche d wind model 
flPudritz &; Norman|[l986{ |Shu et al |["l994[ ), then the re- 
sults imply that these disks are not aligned with the fields 
in the cores from which they formed, or that they are 
preferentially anti-aligned. 

3. THE PROTOSTELLAR DISK FRACTION 

The final input required to calculate the disk fraction is 
a calculation of where disks are expected to form in the pa- 
rameter space of magnetic field strength and field/rotation 
misalignment angle. The most complete numerical study 

5 Telescope Array Doing POLarization 
http: //tadpol . astro . Illinois . edu/ 
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■ Keplerian disk 

▲ Sub-Keplerian disk 

• No disk 



Fig. 1. — The parameter space of ($/M, fi) that describes the 
strength of core magnetic fields, and their orientation relative to 
the angular momentum vectors; /i is the cosine of the misalignment 
angle. Observations indicate that cores are uniformly distributed 
in the r ange <fr/M = 0—1 , and /i = — 1. Points show simulation 
results jjoos et al.||2012| ), indicating whether a simulation with 
those parameters formed no disk, a sub-Keplerian disk, or a 
Keplerian disk. The shaded regions show the range in parameter 
space over which Keplerian disks (red) and non-Keplerian disks 
(blue) form, under our minimal assumptions. The fraction of cores 
that will produce Keplerian disks is the area of the red region, and 
the fraction that will produce any disk at all is the sum of the 
areas of the red and blue regions. 



■ Keplerian disk 
▲ Sub-Keplerian disk 
• No disk 
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Fig. 2. — Same as Figure [l] but now the shaded regions 
show maximal rather than minimal assumptions about where in 
parameter space disks form. 



to date (Joos et al. 2012[ ) includes 18 simulations us- 
ing ideal MHL) that sample values of <&/M in the range 
0.06 — 0.5, and misalignment angles 8 from — 90°. Fig- 
ures [l]and[2] summarize the numerical results, where each 
result is classified as producing a Keplerian disk (^100 
AU in size or larger), a sub-Keplerian disk (again at least 
^100 AU in size), or no disk at all. 

The observed distribution of field strength implies that 
cores are uniformly distributed in <E>/Af from to 1. If 
we adopt the more conservative TADPOL result that 
field/rotation alignment is random, as opposed to pref- 



erentially anti-alignned, then cores are also uniformly 
distributed in fi — cos 9 from to 1, implying that the 
distribution in the ($/M, fi) plane is uniform from to 1 
in both dimensions. In this case the fraction of systems 
that possess disks, Keplerian or otherwise, is simply the 
fraction of the area of the unit square in the (&/M, pt) 
plane over which disks form. To estimate this, we can 
make either minimal or maximal assumptions about disk 
formation based on the simulations. In the minimal 
case we assume that if a disk forms in a simulation at a 
point (^>/M,/j,)i in parameter space, but not at a point 
(<&/M,/i) 2 where either (<J>/M) 2 > ($/A/)i or fi 2 > fJ-i, 
then disks form only at $/M < ($/M)i and /i < /ii. 
A maximal assumption is the opposite: disks form for 
all 3>/M < ($/Af) 2 and < [i 2 . The shaded regions in 
Figures [l] and [2] summarize the two cases. With either 
assumption we can compute the disk fraction simply by 
integrating over the shaded regions. Doing so we find that 
Keplerian disks are expected to form in (11%, 48%) of 
cores, and either Keplerian or non-rotationally-supported 
pseudo-disks in (29%, 50%) of cores, where the first num- 
ber corresponds to the minimal assumption, and the 
second to the maximal. 

4. DISCUSSION AND CONCLUSIONS 

Our calculation that 10 — 50% of cores should produce 
Keplerian disks provides a natural explanation for the 
significant number of Class and I sources around which 
disks have been observed. At present we lack a full census 
of the disk fraction at these early stages, and thus it is 
unclear if the percentages we compute are consistent with 
the actual fraction of such systems that have disks, or if 
some other mechanism will be required. 

While our results help alleviate the problem at the 
Class and I stage, a disk fraction of 50% still too small 
in comparison to what is observed at the Class II stage, 
wher e the disk fraction approaches unity (Haisch et al. 
2001). The question, then, is what mechanism might ex- 
plain such a high disk fraction at these later stages. One 
possibility is that we have been too conservative in adopt- 
ing a uniform distribution of field/rotation orientations, 
and that in fact fields are preferentially anti-aligned with 
core angular momentum vectors. If this is the case, then 
cores will preferentially occupy the region near /i = 0, and 
the disk fraction will be higher than we have estimated. 
This possibility may be checked by fu rther polar i zation 
measurements of the sort performed by Hull et al. (2012 1. 

A second possible explanation for the high traction 
of disks in Class II sources is that it is the result of a 
combination of misalignment with non-ideal MHD effects. 
The values illustrated in Figures [T] and [2] do not include 
any non-ideal MHD effects, such as ion-neutral drift, 
Ohmic dissipation, or turbulent reconnection. Any of 
these effects would probably enhance the ability of disks 
to form, since they would reduce the ability of magnetic 
fields to extract angular momentum from infalling gas. 

A final possibility is that a reduction in the inertia of the 
envelope might yield a Class II disk fraction that exceeds 
that found in Class and Class I sources. In this case the 
disk fraction would only be ^10 — 50% during the Class 
0/1 phase, but would rise to nearly 100% at the transition 
to the Class II phase as the envelope depletes. One may 
distinguish between this scenario and the previous two by 
measuring the disk fraction for Class 0/1 systems. If the 



4 



results are ~10 — 50%, that is consistent with non-ideal 
effects being unimportant and with a uniform distribution 
of disk-rotation angles, in which case the entire difference 
between Class 0/1 and Class II sources would arise from 
a reduction in envelope inertia. On the other hand, a 
disk fraction near unity for Class 0/1 sources would favor 
either non-ideal effects or preferential field/rotation anti- 
alignment as an explanation. 
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